Spatial non-cyclic geometric phase in neutron interferometry 
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We present a split-beam neutron interferometric experiment to test the non-cyclic geometric phase 
tied to the spatial evolution of the system: the subjacent two-dimensional Hilbert space is spanned 
by the two possible paths in the interferometer and the evolution of the state is controlled by phase 
shifters and absorbers. A related experiment was reported previously by some of the authors Q to 
verify the cyclic spatial geometric phase. The interpretation of this experiment, namely to ascribe 
a geometric phase to this particular state evolution, has met severe criticism The extension to 
non-cyclic evolution manifests the correctness of the interpretation of the previous experiment by 
means of an explicit calculation of the non-cyclic geometric phase in terms of paths on the Bloch- 
sphere. The theoretical treatment comprises the cyclic geometric phase as a special case, which is 
confirmed by experiment. 
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I. INTRODUCTION 

Since the discovery of a geometric effect by Berry 3] in 
the shape of an additional phase factor after an adiabatic 
and cyclic transport of a quantum system Berry's phase 
has been intensively investigated and generalized: the ex- 
tension to degenerate subspaces by Wilckzek the re- 
moval of the adiabatic constraint by Aharonov and Anan- 
dan and the cyclic condition by Samuel and Bhandari 
[tII using the early ideas of Pancharatnam f3| and the 
kinematic approach to geomeric phases by Mukunda and 
Simon .Sj, to name a few. In all these contexts the ge- 
ometric phase is dependent only on the geometry of the 
subjacent Hilbert space, but not on the particular dy- 
namics of the system under consideration. Furthermore, 
Manini and Pistolesi j^j proposed an off-diagonal geo- 
metric phase to exhibit the geometry of state space in 
situations where the usual (diagonal) geometric phase is 
undefined. This has been verified experimentally by some 
of the authors (lOj . 

In course of the development of quantum mechanics it 
has become clear that the concept of pure states is not 
sufficient when taking environmental influences causing 
decoherence effects into account. Then one has to use the 
concept of mixed states. Probably the first treatise of a 
geometric phase for mixed states is due to Uhlmann 
in a quantum algebraic context. Another definition of a 
mixed state geometric phase is due to Sjoqvist et al. 
using an interferometric approach for its definition. For 
both definitions one has to keep in mind that there exist 
points in parameter space for which the mixed state ge- 
ometric phases remain undefined provoking an extension 
to off-diagonal mixed state geometric phases [Tsf . 



'Electronic address: sfilipp@ati.ac.at 
tElectronic address: hasegawa@ati.ac. at 



The geometric phase is associated with an evolution of 
a system governed by an Hamiltonian, e. g. a neutron in 
a magnetic field where the geometric phase arises by the 
spinor evolution due to the coupling with the magnetic 
field. Here, we observe a geometric phase as an effect of 
the change in the spatial degrees of freedom in an interfer- 
ometry setup. A proposal to verify the spatial geometric 
phase is due to Sjoqvist using polarized neutrons by 
reversing the roles of the magnetic field and the spatial 
degrees of freedom. Moreover, an experiment using un- 
polarized neutrons has been performed by Hasegawa et 
al. to tsst the cyclic spatial geometric phase by in- 

ducing a relative phase shift of 27r between the interfering 
neutron beams in a perfect silicon single-crystal interfer- 
ometer. The geometric interpretation of this experiment 
has been dismissed by Wagh Q demanding for further 
investigations, namely in the non-cyclic case, which is 
the purpose of the current article. 



II. GEOMETRIC PHASES 

Let us briefly review the basic concepts of geometric 
phases: A geometric phase is a quantity which is deeply 
connected to the curvature of some underlying (state- 
or parameter-) space. A two-dimensional plane in three 
dimensional real space does not have an intrinsic cur- 
vature, but when considering a sphere embedded in eu- 
clidean real space we have to take the curvature of this 
manifold into account. In geometry this curvature is re- 
flected for example in the angle difference of a vector 
transported around a loop along geodesies, i. e. great 
circles: If a vector is pinned onto a sphere and then trans- 
ported along a meridian to the equator, for some angle 
a along the equator and back to the initial point with- 
out changing its length and its direction in the tangent 
plane to the surface of the sphere, the vector will point 
in a different direction with a relative angle of a as the 
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holonomy associated with the loop. If we do the same on 
a two-dimensional plane the initial and the final vector 
will point in the same direction. 

Berry Q was the first who addressed this issue in 
quantum mechanics: He considered a system initially 
in an eigenstate \n{R{t)))t=Q of the governing Hamilto- 
nian H{R{t)) dependent on the parameters R{t) changing 
with time t. As a demonstrative example one may con- 
sider a neutron coupling to a magnetic field H{R{t)) = 
—fl- B{R{t)) due to its magnetic dipole moment /2 = /inf?, 
where a = {(Jx,^y,<^z} are the Pauli matrices and the 
magnetic moment fin of a neutron. Suppose now that 
the neutron is initially polarized in the direction of the 
magnetic field. If the direction of the magnetic field is 
changed adiabatically, i. e. slowly enough to avoid tran- 
sitions to an orthogonal state, the system will stay in the 
eigenstate \n{R{t))) at all times t. Furthermore, when 
tracing out a loop in parameter space the final state 
|\E'(t)) at time r will be the same as the initial state 
up to an additional phase factor: 

\^^^{T)) = e"'''' \n{R{T))) 

The first phase value (pd — ~\ En{R{t))dt is de- 
pendent on the time needed to transverse the loop and 
on the energy En{t) = {t)\H {t)\^ {t)) of the system, 
whereas the second phase (j)g — i §u{n{R)\^ ^\n{R))dR is 
dependent only on the circuit integral in parameter space 
revealing the geometric structure. The latter is termed 
Berry phase or more general geometric phase in contrast 
to the former dynamical phase (j)d. 4>g can be rewritten 
as a surface integral by use of Stoke's Theorem yielding 
(j)g = — Im JpdSVn{R), where F is the surface enclosed 

by the loop in parameter space with dS denoting the area 
element and Vn = V x {n\Vn) in an obvious abbreviated 
notation. For the neutron example - or more generally 
for any spin-1/2 particle - (j)g equals half of the solid angle 
enclosed by the loop as seen from the degeneracy point 
\R\ = in parameter space. This can also be related to 
the example from geometry above where the holonomy 
after the transport of the vector pinned initially to the 
north pole of a sphere equals the solid angle as seen from 
the origin of the sphere. 

Several restrictions have been relaxed in course of the 
years, e. g. extensions to nonadiabatic 0, noncyclic and 
nonunitary 0, and nonpure [Til IT^ geometric phases 
have been made. Important in our case are the gener- 
alizations to the nonadiabatic regime by Aharonov and 
Anandan and to noncyclic paths by Samuel and Bhan- 
dari. 

In this case we have to introduce the Projective Hilbert 
space (Ray space) TZ by identifying all state vectors in 
Hilbert space Ti. which differ only by an overall phase 
factor: 

|0)^|0'):|0'>=e-|0), aeM. (2) 



The stress is therefore shifted from the parameter space 
of the Hamiltonian in case of Berry's construction to state 
space. We are not interested in the changes of the driving 
parameters (as the direction and strength of the magnetic 
field) but in the changes of the state itself. In Berry's con- 
siderations these two spaces are identical since the state 
follows the changes in parameters due to the adiabaticity 
condition. 

In the construction by Aharonov and Anandan one 
considers an open path in Hilbert space which is pro- 
jected to a path in Ray space by use of the equivalence 
relation in Eq. J^J, i. e. the curve C : t G [0,t] i-^ 
|0(i)) e H is projected to C : t & [0, r] i-> 7r(|0(t))) = 
\<j){t)) {(f>{t)\ e 7^ with 10(0)) - |</>(t)). For C an abso- 
lute phase factor of \(l){t)) is immaterial, since the curve 
C is defined via the evolution of the projection operator 
\(j){t)){(j){t)\. The geometric phase is a property of Ray 
space, where C is closed due to the equivalence of the 
initial and the final state |0(O)) ^ |0(t)) and can be cal- 
culated via a surface integral over the area enclosed by 
C. 

One can find many different curves C, C", . . . in Hilbert 
space differing by a phase factor e*"^*-' and yielding the 
same curve in Ray space under the projection map tt. On 
the other hand for a given curve in Ray space there ex- 
ists one distinct curve in Hilbert space fulfilling the par- 
allel transport conditions, namely that two neighbour- 
ing states \4>{t)) and |0(i -I- dt)) in H have the same 
phase, that is to say, {(j){t)\(f){t + dt)) is real and positive. 
This imphes by Taylor expansion that {(f>{t)\-^\(l){t)) = 
[l8j |. For this curve the dynamical phase vanishes as 
one can verify by inserting the Schrodinger equation 
H{t)\<j){t)) = ihd/dt\(j){t)) into the parallel transport con- 
dition. 

The concept can be extended to apply to open paths 
in Ray space where |<^(t)) 7^ |0(O)) by closing the curve 
by a geodesic, i. e. a path in Ray space with the shortest 
distance from |0(r)) ((/)(t)| to \(t){Q)){(j){Q)\. Then one ob- 
tains a well-defined surface area enclosed by the path gen- 
erated by the evolution of the system plus the geodesic 
closure. This surface provides an expression for the ge- 
ometric phase, which has been proven by Samuel and 
Bhandari 0. 

To sum up, for a general evolution of a quantum state 
the state obtains a dynamical phase dependent on the 
energy and time as well as a geometric phase only de- 
pendent on the subjacent geometry of state space. For 
special Hamiltonians which fulfill the parallel transport 
conditions the dynamical phase vanishes, which is also 
the case when the state is transported along a geodesic. 
An example of the latter is an evolution along a great 
circle on a sphere for a two-level system which we will 
encounter in the forthcoming discussion. 
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III. INTERFEROMETRIC SETUP 

Due to Feynmann the description of any two-level 
quantum system is equivalent to the description of a 
spin-1/2 particle. Exploiting this equivalence there is in 
principle no difference between manipulations in the spin 
space of neutrons with the orthogonal basis {| t), | i)} as 
eigenstates of a^-, and momentum space with 
as orthogonal basis vectors corresponding to two direc- 
tions of the neutron beam in an interferometer. In both 
cases one can assign a geometric phase to the particular 
evolution of the initial state. An even more appropriate 
description for the interferometric case for the forthcom- 
ing discussion is in terms of "which-way" basis states 
namely, if the neutron is found in the upper 
beam path after a beamsplitting plate it is said to be in 
the state or in the state Ip"""), if found in the lower 
beam path. In case of a 50 : 50 beamsplitting of the 
incident (neutron) beam into a transmitted beam and a 
reflected beam the associated wave vector after the beam- 
plitter can be written as an equally weighted coherent su- 
perposition of the two paths \q{S)) = l/\/2(|p) +e^^\p^)) 
with the relative phase (5 e R depending on the particular 
physical realization of the beamsplitter. 




FIG. 1: Experimental setup to test the spatial geometric 
phase in a neutron interferometer 

For testing the spatial geometric phase we use a 
double-loop interferometer (cf. Fig. where the in- 
cident unpolarized neutron beam 1-0) is split up into a 
diffracted reference beam \ipref) and a transmitted beam 
\ipt)- The transmitted beam is subjected to further evo- 
lution in the second loop of the interferometer by use 
of beamsplitters (BSl and BS2), an absorber (A) with 
transmission coefficient T and the phase shifter PS2 gen- 
erating a phase shift of 6*"^^ on the upper and e''^^ on the 
lower beam path, respectively, yielding the final state 
IV'/) = Ultpt) = U\p). Here |V't) = \p) since before the 
beam splitter BSl the beam is clearly localized as seen 
from the second loop of the interferometer. The unitary 
matrix U = f/(T, (f>i,(j)2) comprises all the manipulations 
in the second loop: 



IV*) ^ i=(b) -fb^)) A, _L(|p) 



PS2 1 



(e^^Mp) + VTe'^^lp^)) ^ C/|Vt) = |V/{3) 



The geometric phase can then be extracted from the ar- 
gument of the complex valued scalar product between 
the initial and the final state arg('0t|'0/) (when removing 
dynamical contributions as will be discussed later) . This 
is where the reference beam comes into play: IVre/) is 
not subjected to any further evolution, but is stationary 
apart from adding a phase factor e"' by use of the phase- 
shifter PSl. IVre/) propagates towards the beamsplitter 
BS2 from the upper path, thus we can assert it to be in 
the state e*'' \p) . Then by the variable phase shift e*'' one 
can measure the shift of the interference fringes reflecting 
the phase difference between IVre/) and \tpf). 

This preparation of the states is followed by the re- 
combination of the two beams \ipf) and \'ipref) at the 
beamsplitter BS2 and the detection at the detector Do 
in the forward beam. This step can be described by the 
application of the projection operator \q){q\ = l/2(|p) + 
\p'^)){{p\ + {p'\) (with 6 = 0, which can always be achieved 
by an appropriate choice of the phase of the basis states) 
to IV/) as well as to IVre/): 



(4) 



IV}) = \q){qm=K{e^*^+Vfe^^^)\q) 

Wref) - \q){qHref)^K\q), 



where K is some scaling constant. 

The intensity / measured in the detector Dp is pro- 
portional to the absolute square of the superposition 

|V}) + e'''|V;e/): 

\ie^^ + e^f^Vfe^f^)\q)\' = (VWKe/) + i^'f) + 



oc 



-2|(V;e/l^})|cos(r/-arg(V;,/|V})) 



(5) 



We notice a phase shift of the interference pattern by 
arg(V're/IV'/)- This phase shifts corresponds to the Pan- 
charatnam connection between the state |V/) and the 

state \ipref) — \l){l\''Pt) = \q) from which we can extract 
the geometric phase. Explicitly we obtain 



4, = arg(V^.e/|V}; 



arctan 



tan 



T 



,(6) 



where Ac/) 



The geometric phase is defined as 



0g = arg(V,',e/IV'/) -(/-d, (7) 

where (j)d denotes the dynamical part. From Refs. 
and we know that the dynamical part stemming from 
the phase shifter PS2 is given by a weighted sum of the 
phase shifts and 02 with the weights depending on the 
transmission coefficient T. In particular we have 



1 + T 



which vanishes by an appropriate choice of phase shifts 
and transmission, i. e. 0d = for 01 /02 = —T. 

By varying the relative phase A0 from to 27r and 
setting 0^ = the geometric phase 0g can be plotted 
over A0 (cf. Figure O. 
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FIG. 2: Geometric phase (pg [rad] in dependence on the rela- 
tive phase shift [rad] for T = 1/8. 



IV. BLOCH-SPHERE DESCRIPTION 

For every two-level system we can use the Block- sphere 
for depicting the state vectors and evolutions thereof as 
points and curves on a sphere. Then the results obtained 
above, i. e. the shift of the interference pattern in Eq. 

without dynamical contributions, should be equal to 
the (oriented) surface area enclosed by the paths of the 
state vectors on the Bloch-sphere, or, equivalently, to the 
solid angle traced out by the state vectors as seem from 
the origin of the sphere. 

As we can observer in Figure 13 the north pole of the 
sphere can be identified with a state with well known 
path, i. e. an eigenstate of the observable \p){p\- After 
the beam splitter BSl the state \ipt) evolves to an equal 
superposition of upper path and lower path, therefore the 
evolution on the bloch sphere is given by a geodesic from 
the north pole to the equatorial line (the particular point 
on the equator is arbitrary due to the arbitrary choice of 
the phases of the basis vectors). 




FIG. 3: Path of the state in an interferometer on the Bloch 
sphere representing the 2-level system (upper path \p){p\ and 
lower path \p'){p'\). 

The absorber changes the weights of the superposed 
basis states, in particular for the extremal values of T 
parameterized by the angle 9 with T = tan^ 9/2 we end 
up either again with an equally weighted superposition 



for no absorption (T = 1 or = 7r/2) or the state is 
now on the north pole for total absorption (T = or 
9 = 0), since in the latter scenario we know the particle 
has taken the upper path when detecting a neutron in 
Do- For T S (0, 1) the state is encoded as a point on the 
geodesic from the north pole to the equatorial line. 

Due to the phase shifter PS2 we obtain a relative phase 
shift between the superposing states of A(j) = (1)2 — (j^i- 



V2 



i\p) + \P^)) 



V2 



i\p) 



(9) 



This can be depicted as an evolution along a circle of 
latitude on the Bloch sphere with perodicity of 27r. 

The recombination at BS2 followed by the detection of 
the forward beam in Do is represented as a projection to 
the starting point on the equatorial line, i. e. we have to 
close the curve associated with the evolution of the state 
by a geodesic to the point \q){q\ on the sphere as discussed 
for non-cyclic paths in Section ^1 As for the reference 
state li^ref) we note that the phase shift of rj has no 
impact on the position of the state on the Bloch sphere, 
it stays at the north pole. Due to the recombination at 
BS2 and the detection the state is also projected to \q){q\ 
contributing to the forward beam incident to the detector 
Do. 

Th e paths are depicted in Figure ^ in d etail for cyclic 
4(a) as well as non-cyclic evolution |4(b) For a relative 
phase difference greater than 7r/2 w e have to take the di- 
rection of the loops into account. In 4(b) the first loop is 
transversed clockwise, whereas the second loop is trans- 
versed counter-clockwise yielding a positive or negative 
contribution to the geometric phase, respectively. 

The numerical calculation of the surface area F 
enclosed by the path transversed by the neutron is 
straightforward by evaluating the solid angle = 
JpSm9 d{A(f>)d9 via a surface integral and using (j)g = 
—Vt/2. For the cyclic case this integral can be solved 
easily by calculating the segment on the sphere accord- 
ing to Figure 4(a) to obtain (pg = —17/2 = tt{cos9 — 1). 
For the non-cyclic case no analytic expression has been 
found to compare the results with the phase shift of the 
interference fringes appearing in Eq. ((SJ. However, the 
numerical results are equivalent and agree with Figure |21 



V. EXPERIMENTAL RESULTS 

For an experimental test of the spatial geometric 
beam we have used the double-loop perfect-crystal- 
interferometer installed at the S18-beamline at the high- 
flux reactor ILL, Grenoble 0|. A schematic view of the 
setup is shown in Figure ^ Before falling onto the skew- 
symmetric interferometer the incident neutron beam was 
coUimated and monochromatized by the 220-Bragg re- 
flection of a Si perfect crytal monochromator placed in 
the thermal neutron guide H25. The wavelength was 
tuned to give a mean value of Aq = 2.715A. To elim- 
inate the higher harmonics we have used prism-shaped 
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|p)(r 




7r\ 




(a) Cyclic evolution 




(b) Non-cyclic evolution 

FIG. 4: paths on the Blochsphere corresponding to the evo- 
lution of the state in the splitbeam experiment 



silicon wedges. The beam cross-section was confined to 
5 X 5mni^ and an isothermal box enclosed the interfer- 
ometer to achieve reasonable thermal environmental iso- 
lation. For the phase shifters parallel sided aluminium 
plates have been used: 4mm inserted in the first loop 
as PSl and 4mm and 0.5mm, respectively, in the sec- 
ond loop for PS2 yielding a ratio of 1/8 for (t>i/(f)2- To 
avoid dynamical phase contributions a gadoldinium so- 
lution with T — 0.118 [3 has been used as an absorber. 
For a comparison with the theoretically predicted val- 
ues one has to keep in mind that the contrast reflecting 



the coherence properties is different between each of the 
beams in the interferometer. Accounting for this experi- 
mental fact in the theoretical derivation of the geometric 
phase we notice a slighlty flattened curve in Figure |S1 
compared to Figure |21 Nevertheless, one can recognize 
the increase in geometric phase for G [0,7r/2] due to 
the positively oriented surface followed by a decrease due 
to appearance of a counter-clockwise transversed loop on 
the sphere yielding a negative phase contribution. 




FIG. 5: Experimental verification of the spatial geometric 
phase using a neutron interferometry setup 



VI. CONCLUSIONS 

In summary we have shown that one can ascribe a geo- 
metric phase not only to spin evolutions of neutrons, but 
also to evolutions in the spatial degrees of freedom of 
neutrons in an interferometric setup. This equivalence is 
evident from the description of both cases via state vec- 
tors in a two dimensional Hilbert space. However, there 
have been arguments contra the experimental verifica- 
tion in [H which we believe can be settled in favour of a 
geometric phase appearing in the setup described above. 
The twofold calculations of the geometric either in terms 
of a shift in the interference fringes or via surface inte- 
grals in an abstract state space allows for a geometric 
interpretation of the obtained phase shift. 
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